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. ^_, Abstract 

H , An-energy dependent variation of the X-ray emission from the archetypical "Compton- 

thick" Seyfert 2 galaxy NGC 1068 has been observed between two BeppoSAX ob- 
servations, performed one year apart. This variation (at the 2.6(7 level) is highest 
in the 3-10 keV band and declines towards lower and higher energies. The most 
straightforward explanation is a variation of the primary nuclear continuum, which 
is scattered along our line of sight by an electron plasma. If this is indeed the case, 
this evidence allows us to obtain the first direct estimate of the location of the 
scattering medium, at ~1 pc from the nucleus. 
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1 Introduction 



Compton-thick Seyfert 2 galaxies, i.e. sources in which the nucleus is obscured 
by matter with column densities exceeding ~ a~[fl ~ 10^^ cm~^, can be ob- 
served, in X-rays, only in scattered light. In the best known case, NGC 1068, 
at least two scattering components are known to be present (e.g. Iwasawa, 
Fabian & Matt 1997; Matt et al. 1997). One originates from almost neutral 
material, possibly the inner surface of the molecular torus, envisaged by the 
Seyfert unification scenarios, Antonucci & Miller 1985; Antonucci 1993. It is 
recognizable by the characteristic shape of the continuum, peaking between 20 
and 30 keV, which is determined by the competition between photoabsorption 
and Compton scattering. A prominent iron line at 6.4 keV comes along with 
it (George & Fabian 1991; Matt, Perola & Piro 1991). The other originates 
from ionized matter, and is a fainter replica of the nuclear radiation, at least 
as long as self-absorption effects are negligible. Lines from ionized elements, 
most commonly the He- and H-like stages, are superposed to the latter. 

Compton-thick Seyfert 2s are therefore ideally suited for studying the circum- 
nuclear matter in AGN. From X-ray spectroscopy, much can be learned on 
the physical and chemical properties of the reflecting matter (e.g. Netzer 1996; 
Matt, Brandt & Fabian 1996, Netzer & Turner 1998). Even if a full exploita- 
tion of the line diagnostic capabilities should await high spectral resolution 
instruments, first results have already been obtained with moderate resolu- 
tion detectors (Netzer & Turner 1998; Guainazzi et al. 1999, G99 hereinafter). 
To map the location of the reflectors, instead, variability measurements are 
required. While an accurate mapping requires a good sampling performed 
with a high sensitivity detector (and even in this case a real reverberation 
technique cannot be applied due to the fact that the nuclear radiation is not 
directly observable), basic pieces of information can be derived by the crude 
comparison of observations made at different epochs, preferentially with the 
same instrument to avoid cross-calibration problems. 

The Italian-Dutch X-ray satellite BeppoSAX (0.1-200 keV, Boella et al. 1997a) 
observed NGC 1068 twice, about one year apart. We therefore searched for 
variability between the two observations, and found an energy dependent one. 
In Sect. 2 a detailed description of the data reduction is presented, and pos- 
sible systematic effects are discussed in Sect. 3. As the variability survived 
all our checks, in Sect. 4 we discuss the results and show that a possible ex- 
planation of the broadband X-ray variability is in term of a flattening of the 
nuclear component, mirrored by the ionized reflector. In Sect. 5 we derive some 
geometrical and physical properties of the warm scatterer (and of the other 
spectral components) implied by our results. 



Table 1 

NGC 1068 BeppoSAX observation log 



ID 



Start Time 
(UTC) 



ErrdTime T^E^S t^ecs t™/ 

(UTC) (ks) (ks) (ks) 



December 1996 30/12/1996 08:47:25 03/01/1997 05:27:50 61497 100150 62493 
January 1998 11/01/1998 09:52:36 12/01/1998 08:07:50 15408 37331 17657 



2 Observation and data reduction 



NGC 1068 was observed by BeppoSAX on December 1996 and January 1998 
(see Tab. 1). Of the four co- aligned instruments, comprising the BeppoSAX 
scientific payload, we will deal in this paper with data from the Low En- 
ergy Concentrator Spectrometer (LECS, 0.1-10 keV, Parmar et al. 1997), the 
Medium Energy Concentrator Spectrometer (MECS, 1.8-10.5 keV, Boella et 
al 1997b) and the Phoswitch Detector System (PDS, 13-200 keV; Frontera 
et al. 1997). LECS and MECS are imaging gas proportional counters, which 
use identical concentrator systems, with moderate energy resolution (~8% at 
6 keV). The MECS effective area in the two units configuration discussed in 
this paper is about twice as the LECS one in the overlapping energy band. 
The PDS consists of four independent crystal units arranged in pairs, each 
mounted on a rocking collimator to achieve a continuous monitoring of the 
background with a duty-cycle of 96 s. 

Data reduction followed standard procedures, as described in G99. In partic- 
ular. Good Time Intervals (GTI) for scientific products accumulation were 
defined according to the following selection criteria: a) the angle between 
the pointing directions and the Earth's limb was higher than 4°; b) passages 
through the South Atlantic Geomagnetic Anomaly were excluded; a) the an- 
gle between the pointing direction and the direction of the Sun was greater 
than 60°. Count rates were extracted from circular regions of radius 8' and 4' 
around the NGC 1068 centroid in the LECS and MECS, respectively. No other 
point source is detected in these areas at a level higher than 6 x 10"^^ and 
1.6 X 10^^^ erg cm~^ s~^ (3 a local background fluctuations) in the 0.1-2 keV 
and 2-10 keV energy bands, respectively. For comparison, G99 report average 
fluxes of 1.1 and 0.5 x 10~^^ erg cm~^ s~^ for the BeppoSAX observations of 
NGC 1068 in the same bands. PDS data were additionally screened by remov- 
ing flve minutes intervals after any SAGA passage, to allow gain recovery to 
nominal values after instrumental switch-offs. 

We have performed an analysis of several BeppoSAX Crab observations to 
estimate if aging effects in the detectors could produce a secular variation of 
the observed count rate. No effect of this kind is observed, the upper limit 



Table 2 

LEGS background subtracted count rates (in units of 10^^ s^^) when the three dif- 
ferent techniques of background subtraction after Parmar et al. (1999) are adopted 
(details in text) 



0.1-1 keV 1-3 keV 3-10 keV 
December 1996 

BLS 6.29 ±0.10 4.07 ±0.08 1.51 ±0.07 

FOV 6.58 ±0.11 4.02 ±0.09 1.37 ±0.07 

RM 6.28 ±0.11 4.00 ±0.09 1.55 ±0.08 
January 1998 

BLS 5.82 ±0.20 4.01 ±0.17 1.60 ±0.13 

FOV 6.10 ±0.20 3.96 ±0.18 1.45 ±0.13 

RM 5.81 ±0.20 3.94 ±0.18 1.64 ±0.13 

being 0.3% and 0.7% for the LEGS and MEGS, respectively. The count rates 

obtained from the two NGG 1068 observations may be hence directly com- 
pared. 



3 LEGS data reduction issues 



One of the most critical issues in estimating the variability of a relatively faint 
source as NGG 1068 is the background subtraction. In Tab. 2 we compare the 
count rates in three energy bands (0.1-1 keV; 1-3 keV; 3-10 keV) when the 
three background subtraction techniques described in Parmar et al. (1999) are 
employed. They are: a) background extracted from blank sky deep exposures, 
accumulated by the BeppoSAX Science Data Genter (SDG) in the first three 
years of operative life of the mission (BLS); b) background extracted from two 
semi-annuli in the LEGS field of view, suitably renormalized to the expected 
counts under the source extraction region (FOV); c) cosmic background esti- 
mated from the 7-bands ROSAT/PSPG All Sky Survey source-removed count 
rates (RM), and added to the instrumental background. The results differ 
by about 5% , 2% and 12% in the three bands, respectively. This differences 
might be due to local granularity in the contribution of the cosmic background, 
which is globally mapped on scales not lower than 2° (Snowden et al. 1995). On 
the other hand, the instrumental LEGS background is known to exhibit sec- 
ular variations with a dynamical range of about a factor two (T.Oosterbroek, 
private communication). It is therefore impossible to tell which of the three 
techniques is giving the "true" result (see also the discussion in Parmar et 
al. 1999). Although the count rates uncertainties are dominated by system- 



Table 3 

MECS background subtracted count rates (in units of 10^^ s^"*^) 

1.5-3.0 keV 3.0-5.0 keV 5.0-10.0 keV 

December 1996 

2.40 ±0.05 1.44 ±0.04 2.36 ± 0.05 

January 1998 

2.39 ±0.08 1.63 ±0.07 2.57 ± 0.09 

atic effects, the ratio between the count rates observed in the two BeppoSAX 
observations is independent of the background subtraction method used. As- 
suming the BLS one (the same adopted normally in spectral fitting), the ratio 
of the January 1998 versus the December 1996 count rates are: 0.93 ± 0.04, 
and 0.99 ±0.05 and 1.06 ±0.11 in the 0.1-1 keV, 1-3 keV and 3-10 keV energy 
bands, respectively. 

Count rate variations by a few percent could be ascribed to the a different 
positions of the source in the LEGS field of view. It is in facts crossed by 
a mesh of fine grids, which cause a ~5% modulation of the count rates on 
scales ~10", and by the strongback support rib, which causes a ^30% ab- 
sorption at regular distances of 2' (along both the X and Y axis of the focal 
plane). The former effect has a much lower spatial scale than the PSF width 
(~4' at 0.28 keV; ~0.9' at 6 keV), and is therefore averaged out when wide 
enough extraction radii are taken. The latter effect is unlikely to introduce 
a significant variability. The best-fit centroid of the source profiles differ by 
~40" between the two BeppoSAX observations. On this scale, we estimate a 
count rate difference ^ 1.5%, well within the statistical uncertainties of our 
measures. 

Response matrices appropriated for the source position in the two observations 
were extracted and used in the joint LECS/MECS spectral fitting. The results 
of the spectral fittings, described in detail in Appendix below, confirm the 
above outcomes. 



4 MECS data reduction issues 



In Tab. 4 we report the MECS background subtracted count rates for the two 
observations of NGC 1068 in the 1.5-3 keV, 3-5 keV and 5-10 keV energy 
bands. The softest ratio is consistent with unity, while for E > 3 keV the ratio 
is larger than 1 at the 2. 60" level. 

The change in the position of the source within the MECS detector from 



one pointing to the other (constrained by the better MECS statistics to be 
^ 30") imphed a difference in the vignetting ^ 1%. The gain is stable be- 
tween the two observations. The difference between the best-fit centroids of 
the prominent iron hne in the channel space between the two observations is 
0.6 ±0.9 PI (1 channel corresponds to ~50 eV). We note that the MECS does 
not suffer of any support grid mesh obscuration problem. A support ring is 
present at about 10' from the center of the field of view and is almost con- 
centric with it. No change of the MECS effective area on arc minute scales is 
therefore expected, other than that ascribed to the instrumental vignetting. 
Background subtraction has been performed by using blank field observations. 
To verify that the background during the two NGC 1068 observations is stable 
and consistent with that obtained from the blank fields, we have accumulated 
the count rates from the outer region of the detectors, where the contribu- 
tion from the source is negligible. We find that the background count rates 
vary by less than ~ 6% between the two observations. Considering that the 
background count rate contributes ~ 8% to the source plus background count 
rate, we estimate that any spurious source count rate variation induced by 
inaccurate background subtraction will be smaller than ~0.5%. A systematic 
error of 1.5% will be hereinafter conservatively added to the MECS count 
rates to account for possible vignetting estimate and background subtraction 
inaccuracies. 



5 Comparison between the broadband spectra 



In Fig. 1 we show the ratio between the January 1998 and the December 
1996 spectra. This plot includes also the ratio of the PDS count rates in 
the 13-30 keV energy band. As noted before, the count rate ratio is ~0.9 
below 1 keV, becomes consistent with unity in the 1-3 keV band, reaches a 
maximum (1.098 ± 0.038) between 3 and 10 keV and then decreases again, 
returning consistent with 1 above 15 keV. 

The maximum of the variability occurs in the energy range where the contri- 
bution of the warm scattered component is the highest (G99). The amount of 
the observed variability is admittedly small, and the spectral complexity ham- 
pers a clear identification of which component(s) is responsible for it. However, 
a comparison of Fig. 1 with G99 suggests that the observed variability might 
be due to a flattening of the warm scattered component around a pivot point 
located somewhere in the 1-4 keV band. More complex scenarios cannot be 
ruled out, and the interplay of spectral steepness and normalizations can fur- 
ther contribute to confuse the picture. Variability measurements of the lines 
associated with one or both the reflection components should help enlighten- 
ing on the real behavior, but for this purpose much higher spectral resolution 
detectors are required. 



Fig. 1. Count rate ratio between the January 1998 and December 1996 observations 
as a function of energy. Empty squares: LEGS; filled circles: MEGS; cross: PDS. 
The dotted hnes indicate the MEGS measures in the 3-5 keV and 5-10 keV sepa- 
rately. The solid and dashed line in the PDS data point represent the statistical and 
statistical plus systematic uncertainties after Guainazzi & Matteuzzi (1997). The 
3-10 keV LEGS measure (1.06 it 0.11) is consistent with the MEGS one and not 
shown for sake of clarity 
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In Appendix, we show the comparison between the fitting results for the De- 
cember 1996 and January 1998 spectra. 



6 Discussion 



We report on the observation of an energy dependent flux variability in NGC 1068 
on timescales of about one year. The flux ratio of the 1998 over the 1996 ob- 
servations is higher than 1 between 3 and 10 keV, and decreases towards lower 
and higher energies. Although the result is signiflcant at the 2.6a level only, its 
possible implications are rather far-reaching, and we cannot refrain to discuss 
them in this Section. The next generation of high throughput, high resolution 
X-ray detectors will allow a robust testing of the above results and of the 
forthcoming hypothesis. 



General arguments on the spectral deconvolution suggest that the observed 
variation is best explained in terms of a change of the warm reflection com- 
ponent, echoing the same variation in the nuclear spectrum. The statistics 



is not good enough, however, to tell whether it is the spectral index or the 
normalization of the warm-scattered component (or both), that is responsible 
for the observed variability (see Appendix). Similarly, our data cannot test in- 
dependently the hypothesis that self- absorbing effects in the warm-scattered 
component are negligible. 

Apart from the details, it is rather plausible that the measured variation is 
mainly due to the warm reflection component. This may permit to derive 
information on the size and location of the reflecting matter. Light crossing 
arguments suggest that the size of the region cannot exceed one light-year, i.e 
one-third of a parsec. As the light-curve of the nuclear emission is unknown, 
we cannot determine with certainty the distance of the reflector to the nucleus. 
However, an order-of-magnitude estimate may be derived by recalling that 
the 2-10 keV luminosity of the warm reflection component is about 1.5x10'^^ 
erg s^^, and that its Thomson optical depth, tt, is estimated to be a fewx 10^^ 
(G99). The relation between nuclear and reflected luminosity may be written 
as: 



>re/ 



rrfLr. 



where / is a geometrical factor basically equal to the solid angle subtended by 
the reflecting matter to the nucleus, in units of Air. Therefore, if the nuclear 
luminosity is not much larger than 10^^ erg s~^ (Iwasawa et al. 1997), / cannot 
be much lower than 1. If this is the case, the distance of the reflecting region 
from the nucleus cannot be much larger than its size, and should therefore be 
of the order of a parsec or less. This crude estimate puts the warm scattering 
medium in an intermediate location between the Broad (Peterson et al. 1991; 
Stirpe et al. 1994; Dietrich et al. 1998) and the Narrow Line Regions (Axon et 
al. 1997). This supports the idea that this medium coincides with that respon- 
sible for the almost energy-independent scattering of the nuclear continuum, 
observed as the appearance of broad lines in the polarized light optical spectra 
of NGC 1068 (Antonucci & Miller 1985). Because in this source the torus is 
probably observed almost edge-on (Greenhill et al. 1997; Matt et al. 1997), 
the warm reflector, in order to be visible, must be placed outside the torus 
itself. This implies a rather compact torus, with inner side of order of 1 pc 
or less. This is in good agreement with the radio continuum measurements, 
which suggest the presence of a sub-parsec edge-on disk around the NGC 1068 
nucleus (Gallimore et al. 1998). 

Netzer & Turner (1998) and G99 suggested that the bulk of the ionized emis- 
sion lines observed in the soft X-ray spectrum of NGC 1068 originates in the 
same medium. The fact that transitions from He-like stages of elements from 
oxygen to iron are present, suggests that the ionization and/or geometrical 
structure of the scattering medium must be complex. This would be confirmed 
by the fact that, while most of the lines do not change significantly between the 



two observations, the Kq, Nexi varies by a factor two (see Appendix). Given 
the instrumental and statistical limits of our measurements, we refrain from 
pursuing this point any further. This results, however, indicates how profitable 
the study of the line dynamics could be, when the new generation of high res- 
olution/high throughput detectors will be operative. New constraints on the 
physical conditions of the matter in the circumnuclear AGN environment can 
be set if monitoring with a sufficiently fine sampling will be available. 

Comparatively little can be said on the variability of the other spectral com- 
ponents (see Appendix). The Compton-refiected continuum is detected only in 
the longer December 1996 observation. However, the January 1998 upper limit 
on its 5-200 keV fiux is not inconsistent with the December 1996 one. The 
thermal component parameters are remarkably constant. Again, this strength- 
ens the idea that this component is connected to the strong starburst ring of 
approximately 1 kpc size [i.e.: about 15" at the distance of 14.4 Mpc (Tully 
1988)], which protrudes a bar towards the nucleus (Scoville 1988). Actually, 
half of the emission at 0.8 keV comes from an extended region of ~13 kpc 
scale (Wilson et al. 1992). It is therefore straightforward to think that the 
BeppoSAX soft X-ray spectrum is dominated by the same component. 



APPENDIX — Time-resolved spectroscopy: an exercise 



We have tried to separately fit the best model of G99 to the December 1996 and 
January 1998 LECS/MECS/PDS spectra. In principle this procedure could al- 
low us to determine exactly which spectral component (s) is responsible for the 
observed variability. In practise this is a very hard task, although BeppoSAX 
at least has the advantage, over previous missions, of the broad energy cov- 
erage which allows the best modeling of the continuum so far. We report the 
relevant results hereinafter mainly for sake of completeness. We extracted the 
MEGS spectra from a region of 6' to achieve a better signal-to-noise ratio. The 
continuum model is characterized by the spectral index and normalization of 
the warm scattered (F^g, N^s) and Gompton-refiected (r^, Ncr) primary con- 
tinuum components, and by the temperature, abundance and normalization 
of the thermal component (kT, Z, Nth; the mekal model in XSPEC has been 
used). The Galactic contribution to the photoelectric absorption has been held 
fixed to the value 2.9 x lO^'' cm~^, following Murphy et al. (1996), which is con- 
sistent with the value obtained in both the December 1996 and January 1998 
fits if left as a free parameter. The results are summarized in Tab. 4 and Tab. 5 
for the continuum and emission lines, respectively. Uncertainties are shown 
at 1 (T level for two interesting parameters (A^^ = 2.30) for all the spectral 
components, except for the thermal optically thin plasma, for which the values 
corresponding to three interesting parameters are reported (A^^ = 3.53). In 
the December 1996 observation, r„s and r^ turned out to be mutually consis- 



Table 4 

Best-fit results with the baseline G99 model for the December 1996 and January 

1998 BeppoSAX NGC 1068 observations 



Parameter 


1996 December 


January 


1998 


r 


2.09 ±0.16 


i-26±S:i 


1.8 ±0.3 


N„s" (10-^ photons s~i cm~2) 


8.3±i:? 


5.o±l:^ 


7±| 


Tcr 


— -■- ws 


2.lVl.25^ 


2.09'' 


p5-200 kev (10-11 erg s-1 cm-2) 


3.2±{:I 


< 1.7/< 2.9 


3.2^ 


kT (eV) 


390 ± 40 


390±i[j 


380±i[j 


Z(%) 


i.6±C:i 


< 1.3 


< 1.9 


Nth^ (10-2 s-5) 


5.4 ±0.8 


5-7±li 


5-3±i:i 


XV dof 


365.6/333 


150.0/156 


156.2/156 



"normalization at 1 keV of the warm scattered component 

^fixed 

'^normalization of the thermal component 

tent, and therefore we report in Tab. 5 the results when these parameters are 
tied together. The only significant difference between the two fits is that the 
warm scattered component is much flatter and brighter in January 1998, the 
2-10 keV flux increasing from (1.84±[j:f^) to (3.4±?;J) x 10"^^ g^g cm-^ s^^ 
This suggest that it has pivoted around about 1 keV between the two obser- 
vations. 

The Compton reflected component is actually detected only in the December 
1996 observation, which has by far the higher signal-to-noise ratio, even if the 
upper limit of the 5-200 keV flux in the second observation is consistent with 
the value in the flrst one (Table 5). The clear presence of the 6.4 keV iron 
line also in the second observation (see below) strongly argues against the 
disappearing of the cold reflection. If the CR component is forced to be the 
same, in shape and flux, as in the flrst observation, it is actually impossible 
to tell whether it was the spectral index or the normalization of the WS 
component that changed (see column 4 of tab. 5). 

In Tab. 4, we show the best-flt results for the emission lines, when the con- 
tinuum models of columns 2 and 3 of Tab. 4 are employed. Most of them 
do not exhibit a variation of either the best-flt centroid energy or the inten- 
sity/Equivalent Width (EW). The only remarkable exception if the Nexi one, 
whose flux doubles when the warm scattered continuum flattens, yielding an 
increase of the EW by a comparable factor. On the other hand, the ionized 
iron lines do not show any signiflcant variation of the intensity, while the 
underlying warm scattered continuum increased by a factor about two. This 
would imply a decrease of their EW by the same factor, rather surprising if 
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Table 5 

Emission line best-fit parameters. The continuum models are those in columns 2 
and 3 (for Fcr = 1.25) in Tab. 4 for the December 1996 and January 1998 observa- 
tions, respectively. Equivalent Width (EW) are calculated against the warm scat- 
tered continuum only, except the Fei, which is given against the Compton-reflected 
continuum only 



"^^ fixed 



Line 


E 


I 


EW 




(keV) 


(photons cm^^ s^-*^) 


(eV) 


December 1996 






Neix 


0.95±0:gi 


(4.2±}i) X 10"^ 


570±150 


Sixiii 


1.88±0:gi 


(8.6 ± 1.3) X 10-5 


460 ± 70 


Sxv 


2.46 ± 0.04 


(3.1±[]:f) X 10-5 


280±50 


Fei 


6.4'' 


(5.3±[]i) X 10-5 


im±'Z 


Fexxv 


G.T" 


{5.3±\i) X 10-5 


3500 ± 800 


Fexxvi 


6.96" 


(2.1±j]:9) X 10-5 


1500±600 


January 


1998 






Neix 


o.9i±g:gt 


(8±i) X 10-4 


1500±600 


Sixiii 


1.84 ±0.05 


(8±i) X 10-5 


380±i40 


Sxv 


o A c,_|_0.06 


(5.o±}:t) X 10-5 


330 ± 90 


Fei 


6.4'* 


(6.o±}:t) X 10-5 


>2200 


Fexxv 


6.7" 


(5±i) X 10-5 


1400±™[1 


Fexxvi 


6.96" 


(3.4±}:5) X 10-5 


1000±|00 



the continuum and the lines are produced in one and the same medium. No 
problems of this kind are present with the fit of column 4 of Tab. 4, however, 
which is therefore again to be preferred on plausibility grounds. 
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